Genome-Wide Analysis and Expression Profiling of the Glutathione Peroxidase-like Enzyme Gene Family in Solanum tuberosum.

Glutathione peroxidase-like enzyme is an important enzymatic antioxidant in plants. It is involved in scavenging reactive oxygen species, which can effectively prevent oxidative damage and improve resistance. GPXL has been studied in many plants but has not been reported in potatoes, the world's fourth-largest food crop. This study identified eight StGPXL genes in potatoes for the first time through genome-wide bioinformatics analysis and further studied the expression patterns of these genes using qRT-PCR. The results showed that the expression of StGPXL1 was significantly upregulated under high-temperature stress, indicating its involvement in potato defense against high-temperature stress, while the expression levels of StGPXL4 and StGPXL5 were significantly downregulated. The expression of StGPXL1, StGPXL2, StGPXL3, and StGPXL6 was significantly upregulated under drought stress, indicating their involvement in potato defense against drought stress. After MeJA hormone treatment, the expression level of StGPXL6 was significantly upregulated, indicating its involvement in the chemical defense mechanism of potatoes. The expression of all StGPXL genes is inhibited under biotic stress, which indicates that GPXL is a multifunctional gene family, which may endow plants with resistance to various stresses. This study will help deepen the understanding of the function of the potato GPXL gene family, provide comprehensive information for the further analysis of the molecular function of the potato GPXL gene family as well as a theoretical basis for potato molecular breeding.

Allelic complexity of KMT2A partial tandem duplications in acute myeloid leukemia and myelodysplastic syndromes.

KMT2A partial tandem duplication (KMT2A-PTD) is an adverse risk factor in acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS), a potential therapeutic target, and an attractive marker of measurable residual disease. High initial KMT2A-PTD RNA levels have been linked to poor prognosis, but mechanisms regulating KMT2A-PTD expression are not well understood. Although KMT2A-PTD has been reported to affect only a single allele, it has been theorized but not proven that genomic gains of a monoallelic KMT2A-PTD may occur, thereby potentially driving high expression and disease progression. In this study, we identified 94 patients with KMT2A-PTDs using targeted DNA next-generation sequencing (NGS) and found that 16% (15/94) had complex secondary events, including copy-neutral loss of heterozygosity and selective gain involving the KMT2A-PTD allele. High copy numbers indicating complexity were significantly enriched in AML vs MDS and correlated with higher RNA expression. Moreover, in serial samples, complexity was associated with relapse and secondary transformation. Taken together, we provide approaches to integrate quantitative and allelic assessment of KMT2A-PTDs into targeted DNA NGS and demonstrate that secondary genetic events occur in KMT2A-PTD by multiple mechanisms that may be linked to myeloid disease progression by driving increased expression from the affected allele.

A novel TP53 tandem duplication in a child with Li-Fraumeni syndrome.

Li-Fraumeni syndrome (LFS) is one of the most common cancer predisposition syndromes that affects both children and adults. Individuals with LFS are at an increased risk of developing various types of cancer over their lifetime including soft tissue sarcomas, osteosarcomas, breast cancer, leukemia, brain tumors, and adrenocortical carcinoma. Heterozygous germline pathogenic variants in the tumor suppressor gene TP53 are the known causal genetic defect for LFS. Single-nucleotide variants (SNVs) including missense substitutions that occur in the highly conserved DNA binding domain of the protein are the most common alterations, followed by nonsense and splice site variants. Gross copy-number changes in TP53 are rare and account for <1% of all variants. Using next-generation sequencing (NGS) panels, we identified a paternally inherited germline intragenic duplication of TP53 in a child with metastatic osteosarcoma who later developed acute myeloid leukemia (AML). Transcriptome sequencing (RNA-seq) demonstrated the duplication was tandem, encompassing exons 2-6 and 28 nt of the untranslated region (UTR) upstream of the start codon in exon 2. The inclusion of the 28 nt is expected to result in a frameshift with a stop codon 18 codons downstream from the exon 6, leading to a loss-of-function allele. This case highlights the significance of simultaneous identification of both significant copy-number variants as well as SNVs/indels using NGS panels.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors.

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

The impact of whole genome duplications on the human gene regulatory networks.

This work studies the effects of the two rounds of Whole Genome Duplication (WGD) at the origin of the vertebrate lineage on the architecture of the human gene regulatory networks. We integrate information on transcriptional regulation, miRNA regulation, and protein-protein interactions to comparatively analyse the role of WGD and Small Scale Duplications (SSD) in the structural properties of the resulting multilayer network. We show that complex network motifs, such as combinations of feed-forward loops and bifan arrays, deriving from WGD events are specifically enriched in the network. Pairs of WGD-derived proteins display a strong tendency to interact both with each other and with common partners and WGD-derived transcription factors play a prominent role in the retention of a strong regulatory redundancy. Combinatorial regulation and synergy between different regulatory layers are in general enhanced by duplication events, but the two types of duplications contribute in different ways. Overall, our findings suggest that the two WGD events played a substantial role in increasing the multi-layer complexity of the vertebrate regulatory network by enhancing its combinatorial organization, with potential consequences on its overall robustness and ability to perform high-level functions like signal integration and noise control. Lastly, we discuss in detail the RAR/RXR pathway as an illustrative example of the evolutionary impact of WGD duplications in human.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L.

BACKGROUNDS: Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae. RESULTS: In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists. CONCLUSIONS: The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.

In-silico genome wide analysis of Mitogen activated protein kinase kinase kinase gene family in C. sinensis.

Mitogen activated protein kinase kinase kinase (MAPKKK) form the upstream component of MAPK cascade. It is well characterized in several plants such as Arabidopsis and rice however the knowledge about MAPKKKs in tea plant is largely unknown. In the present study, MAPKKK genes of tea were obtained through a genome wide search using Arabidopsis thaliana as the reference genome. Among 59 candidate MAPKKK genes in tea, 17 genes were MEKK-like, 31 genes were Raf-like and 11 genes were ZIK- like. Additionally, phylogenetic relationships were established along with structural analysis, which includes gene structure, its location as well as conserved motifs, cis-acting regulatory elements and functional domain signatures that were systematically examined. Also, on the basis of one orthologous gene found between tea and Arabidopsis, functional interaction was carried out in C. sinensis based on an Arabidopsis association model. The expressional profiles indicated major involvement of MAPKKK genes from tea in response to various abiotic stress factors. Taken together, this study provides the targets for additional inclusive identification, functional study, and provides comprehensive knowledge for a better understanding of the MAPKKK cascade regulatory network in C. sinensis.

Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.).

Cystatins, as reversible inhibitors of papain-like and legumain proteases, have been identified in several plant species. Although the cystatin family plays crucial roles in plant development and defense responses to various stresses, this family in wheat (Triticum aestivum L.) is still poorly understood. In this study, 55 wheat cystatins (TaCystatins) were identified. All TaCystatins were divided into three groups and both the conserved gene structures and peptide motifs were relatively conserved within each group. Homoeolog analysis suggested that both homoeolog retention percentage and gene duplications contributed to the abundance of the TaCystatin family. Analysis of duplication events confirmed that segmental duplications played an important role in the duplication patterns. The results of codon usage pattern analysis showed that TaCystatins had evident codon usage bias, which was mainly affected by mutation pressure. TaCystatins may be regulated by cis-acting elements, especially abscisic acid and methyl jasmonate responsive elements. In addition, the expression of all selected TaCystatins was significantly changed following viral infection and cold stress, suggesting potential roles in response to biotic and abiotic challenges. Overall, our work provides new insights into TaCystatins during wheat evolution and will help further research to decipher the roles of TaCystatins under diverse stress conditions.

Molecular and clinical insights into complex genomic rearrangements related to MECP2 duplication syndrome.

MECP2 duplication syndrome (MDS) is caused by copy number variation (CNV) spanning the MECP2 gene at Xq28 and is a major cause of intellectual disability (ID) in males. Herein, we describe two unrelated males harboring non-recurrent complex Xq28 rearrangements associated with MDS. Copy number gains were initially detected by quantitative real-time polymerase chain reaction and further delineated by high-resolution array comparative genomic hybridization, familial segregation, expression analysis and X-chromosome inactivation (XCI) evaluation in a carrier mother. SNVs within the rearrangements and/or fluorescent in situ hybridization (FISH) were used to assess the parental origin of the rearrangements. Patient 1 exhibited an intrachromosomal rearrangement, whose structure is consistent with a triplicated segment presumably embedded in an inverted orientation between two duplicated sequences (DUP-TRP/INV-DUP). The rearrangement was inherited from the carrier mother, who exhibits extreme XCI skewing and subtle psychiatric symptoms. Patient 2 presented a de novo (X;Y) unbalanced translocation resulting in duplication of Xq28 and deletion of Yp, originated in the paternal gametogenesis. Neurodevelopmental trajectory and non-neurological symptoms were consistent with previous reports, with the exception of cerebellar vermis hypoplasia in patient 2. Although both patients share the core MDS phenotype, patient 1 showed MECP2 transcript levels in blood similar to controls. Understanding the molecular mechanisms related to MDS is essential for designing targeted therapeutic strategies.

Lineage-Specific Genes and Family Expansions in Dictyostelid Genomes Display Expression Bias and Evolutionary Diversification during Development.

Gene duplications generate new genes that can contribute to expression changes and the evolution of new functions. Genomes often consist of gene families that undergo expansions, some of which occur in specific lineages that reflect recent adaptive diversification. In this study, lineage-specific genes and gene family expansions were studied across five dictyostelid species to determine when and how they are expressed during multicellular development. Lineage-specific genes were found to be enriched among genes with biased expression (predominant expression in one developmental stage) in each species and at most developmental time points, suggesting independent functional innovations of new genes throughout the phylogeny. Biased duplicate genes had greater expression divergence than their orthologs and paralogs, consistent with subfunctionalization or neofunctionalization. Lineage-specific expansions in particular had biased genes with both molecular signals of positive selection and high expression, suggesting adaptive genetic and transcriptional diversification following duplication. Our results present insights into the potential contributions of lineage-specific genes and families in generating species-specific phenotypes during multicellular development in dictyostelids.

Genome-wide analysis of the NF-Y gene family and their roles in relation to fruit development in Tartary buckwheat (Fagopyrum tataricum).

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor playing crucial roles in various biological process in plant. However, thorough research on NF-Y gene family of Tartary buckwheat (Fagopyrum tataricum) is little. In this study, 38 FtNF-Y genes (12 FtNF-YAs, 17 FtNF-YBs, and 9 FtNF-YCs) were identified and renamed on the basis of their subfamily and chromosomal location. Their gene structure, genomic mapping, motif composition, conserved domain, phylogenetic relationships, cis-acting elements and gene expression were investigated. Illustration of gene structures and conserved domains of FtNF-Ys revealed their functional conservation and specificity. Construction of phylogenetic trees of NF-Ys in Tartary buckwheat, Arabidopsis, tomato, rice and banana, allowed us to predict functional similarities among NF-Ys from different species. Gene expression analysis displayed that twenty-four FtNF-Ys were expressed in all the tissues and the transcript levels of them were different, suggesting their function varieties. Moreover, expression profiles of twenty FtNF-Ys along five different fruit development stages acquired by real-time quantitative PCR (RT-qPCR) demonstrated distinct abundance diversity at different stages, providing some clues of potential fruit development regulators. Our study could provide helpful reference information for further function characterization of FtNF-Ys and for the fruit quality enhancement of Tartary buckwheat.

Deciphering complex genome rearrangements in C. elegans using short-read whole genome sequencing.

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Evolutionary and functional relationships in the ribosome biogenesis SBDS and EFL1 protein families.

Nascent ribosomal 60S subunits undergo the last maturation steps in the cytoplasm. The last one involves removing the anti-association factor eIF6 from the 60S ribosomal surface by the joint action of the Elongation Factor-like 1 (EFL1) GTPase and the SBDS protein. Herein, we studied the evolutionary relationship of the EFL1 and EF-2 protein families and the functional conservation within EFL1 orthologues. Phylogenetic analysis demonstrated that the EFL1 proteins are exclusive of eukaryotes and share an evolutionary origin with the EF-2 and EF-G protein families. EFL1 proteins originated by gene duplication from the EF-2 proteins and specialized in ribosome maturation while the latter retained their function in translation. Some organisms have more than one EFL1 protein resulting from alternative splicing, while others are encoded in different genes originated by gene duplication. However, the function of these alternative EFL1 proteins is still unknown. We performed GTPase activity and complementation assays to study the functional conservation of EFL1 homologs alone and together with their SBDS counterparts. None of the orthologues or cross-species combinations could replace the function of the corresponding yeast EFL1•SBDS binomial. The complementation of SBDS interspecies chimeras indicates that domain 2 is vital for its function together with EFL1 and the 60S subunit. The results suggest a functional species-specificity and possible co-evolution between EFL1, SBDS, and the 60S ribosomal subunit. These findings set the basis for further studies directed to understand the molecular evolution of these proteins and their impact on ribosome biogenesis and disease.

Identification, characterization and expression analysis of lineage-specific genes within mangrove species Aegiceras corniculatum.

Lineage-specific genes (LSGs) are the genes that have no recognizable homology to any sequences in other species, which are important drivers for the generation of new functions, phenotypic changes, and facilitating species adaptation to environment. Aegiceras corniculatum is one of major mangrove plant species adapted to waterlogging and saline conditions, and the exploration of aegiceras-specific genes (ASGs) is important to reveal its adaptation to the harsh environment. Here, we performed a systematic analysis on ASGs, focusing on their sequence characterization, origination and expression patterns. Our results reveal that there are 4823 ASGs in the genome, approximately 11.84% of all protein-coding genes. High proportion (45.78%) of ASGs originate from gene duplication, and the time of gene duplication of ASGs is consistent with the timing of two genome-wide replication (WGD) events that occurred in A. corniculatum, and also coincides with a short period of global warming during the Paleocene-Eocene Maximum (PETM, 55.5 million years ago). Gene structure analysis showed that ASGs have shorter protein lengths, fewer exons, and higher isoelectric point. Expression patterns analysis showed that ASGs had low levels of expression and more tissue-specific expression. Weighted gene co-expression network analysis (WGCNA) revealed that 86 ASGs co-expressed gene modules were primarily involved in pathways related to adversity stress, including plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, peroxisome and pentose phosphate pathway. This study provides a comprehensive analysis of the characteristics and potential functions of ASGs and identifies key candidate genes, which will contribute to the subsequent further investigation of the adaptation of A. corniculatum to intertidal coastal wetland habitats.

Performances of NIPT for copy number variations at different sequencing depths using the semiconductor sequencing platform.

OBJECTIVE: To evaluate the performance of noninvasive prenatal testing (NIPT) and NIPT-PLUS for the detection of genome-wide microdeletion and microduplication syndromes (MMSs) at different sequencing depths. The NIPT sequencing depth was 0.15X, and the data volume was 3 million reads; the NIPT-PLUS sequencing depth was 0.4X, and the data volume was 8 million reads. METHODS: A cohort of 50,679 pregnancies was recruited. A total of 42,969 patients opted for NIPT, and 7710 patients opted for NIPT-PLUS. All high-risk cases were advised to undergo invasive prenatal diagnosis and were followed up. RESULTS: A total of 373 cases had a high risk of a copy number variation (CNV) as predicted by NIPT and NIPT-PLUS: NIPT predicted 250 high-risk CNVs and NIPT-PLUS predicted 123. NIPT-PLUS increased the detection rate by 1.02% (0.58% vs 1.60%, p < 0.001). A total of 291 cases accepted noninvasive prenatal diagnosis, with 197 cases of NIPT and 94 cases of NIPT-PLUS. The PPV of CNV > 10 Mb for NIPT-PLUS was significantly higher than that for NIPT (p = 0.02). The total PPV of NIPT-PLUS was 12.56% higher than that of NIPT (43.61% vs 30.96%, p = 0.03). CONCLUSION: NIPT-PLUS had a better performance in detecting CNVs in terms of the total detection rate and total PPV. However, great care must be taken in presenting results and providing appropriate counseling to patients when deeper sequencing is performed in clinical practice.

Genome-wide identification of PTI1 family in Setaria italica and salinity-responsive functional analysis of SiPTI1-5.

BACKGROUND: PTI1 (Pto-interacting 1) protein kinase belongs to the receptor-like cytoplasmic kinase (RLCK) group of receptor-like protein kinases (RLK), but lack extracellular and transmembrane domains. PTI1 was first identified in tomato (Solanum lycopersicum) and named SlPTI1, which has been reported to interact with bacterial effector Pto, a serine/threonine protein kinase involved in plant resistance to bacterial disease. Briefly, the host PTI1 specifically recognizes and interacts with the bacterial effector AvrPto, which triggers hypersensitive cell death to inhibit the pathogen growth in the local infection site. Previous studies have demonstrated that PTI1 is associated with oxidative stress and hypersensitivity. RESULTS: We identified 12 putative PTI1 genes from the genome of foxtail millet (Setaria italica) in this study. Gene replication analysis indicated that both segmental replication events played an important role in the expansion of PTI1 gene family in foxtail millet. The PTI1 family members of model plants, i.e. S. italica, Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), maize (Zea mays), S. lycopersicum, and soybean (Glycine max), were classified into six major categories according to the phylogenetic analysis, among which the PTI1 family members in foxtail millet showed higher degree of homology with those of rice and maize. The analysis of a complete set of SiPTI1 genes/proteins including classification, chromosomal location, orthologous relationships and duplication. The tissue expression characteristics revealed that SiPTI1 genes are mainly expressed in stems and leaves. Experimental qRT-PCR results demonstrated that 12 SiPTI1 genes were induced by multiple stresses. Subcellular localization visualized that all of foxtail millet SiPTI1s were localized to the plasma membrane. Additionally, heterologous expression of SiPTI1-5 in yeast and E. coli enhanced their tolerance to salt stress. CONCLUSIONS: Our results contribute to a more comprehensive understanding of the roles of PTI1 protein kinases and will be useful in prioritizing particular PTI1 for future functional validation studies in foxtail millet.

Genomic analyses of new genes and their phenotypic effects reveal rapid evolution of essential functions in Drosophila development.

It is a conventionally held dogma that the genetic basis underlying development is conserved in a long evolutionary time scale. Ample experiments based on mutational, biochemical, functional, and complementary knockdown/knockout approaches have revealed the unexpectedly important role of recently evolved new genes in the development of Drosophila. The recent progress in the genome-wide experimental testing of gene effects and improvements in the computational identification of new genes (< 40 million years ago, Mya) open the door to investigate the evolution of gene essentiality with a phylogenetically high resolution. These advancements also raised interesting issues in techniques and concepts related to phenotypic effect analyses of genes, particularly of those that recently originated. Here we reported our analyses of these issues, including reproducibility and efficiency of knockdown experiment and difference between RNAi libraries in the knockdown efficiency and testing of phenotypic effects. We further analyzed a large data from knockdowns of 11,354 genes (~75% of the Drosophila melanogaster total genes), including 702 new genes (~66% of the species total new genes that aged < 40 Mya), revealing a similarly high proportion (~32.2%) of essential genes that originated in various Sophophora subgenus lineages and distant ancestors beyond the Drosophila genus. The transcriptional compensation effect from CRISPR knockout were detected for highly similar duplicate copies. Knockout of a few young genes detected analogous essentiality in various functions in development. Taken together, our experimental and computational analyses provide valuable data for detection of phenotypic effects of genes in general and further strong evidence for the concept that new genes in Drosophila quickly evolved essential functions in viability during development.

The clinical importance of tandem exon duplication-derived substitutions.

Most coding genes in the human genome are annotated with multiple alternative transcripts. However, clear evidence for the functional relevance of the protein isoforms produced by these alternative transcripts is often hard to find. Alternative isoforms generated from tandem exon duplication-derived substitutions are an exception. These splice events are rare, but have important functional consequences. Here, we have catalogued the 236 tandem exon duplication-derived substitutions annotated in the GENCODE human reference set. We find that more than 90% of the events have a last common ancestor in teleost fish, so are at least 425 million years old, and twenty-one can be traced back to the Bilateria clade. Alternative isoforms generated from tandem exon duplication-derived substitutions also have significantly more clinical impact than other alternative isoforms. Tandem exon duplication-derived substitutions have >25 times as many pathogenic and likely pathogenic mutations as other alternative events. Tandem exon duplication-derived substitutions appear to have vital functional roles in the cell and may have played a prominent part in metazoan evolution.

Protein kinase C family evolution in jawed vertebrates.

Protein kinase C (PKC) was one of the first kinases identified in human cells. It is now known to constitute a family of kinases that respond to diacylglycerol, phosphatidylserine and for some family members, Ca2+. They have a plethora of different functions, such as cell cycle regulation, immune response and memory formation. In mammals, 12 PKC family members have been described, usually divided into 4 different subfamilies. We present here a comprehensive evolutionary analysis of the PKC genes in jawed vertebrates with special focus on the impact of the two tetraploidizations (1R and 2R) before the radiation of jawed vertebrates and the teleost tetraploidization (3R), as illuminated by synteny and paralogon analysis including many neighboring gene families. We conclude that the vertebrate predecessor had five PKC genes, as tunicates and lancelets still do, and that the PKC family should therefore ideally be organized into five subfamilies. The 1R and 2R events led to a total of 12 genes distributed among these five subfamilies. All 12 genes are still present in some of the major lineages of jawed vertebrates, including mammals, whereas birds and cartilaginous fishes have lost one member. The 3R event added another nine genes in teleosts, bringing the total to 21 genes. The zebrafish, a common experimental model animal, has retained 19. We have found no independent gene duplications. Thus, the genome doublings completely account for the complexity of this gene family in jawed vertebrates and have thereby had a huge impact on their evolution.

Reconstructing the evolutionary history of the oxytocin and vasotocin receptor gene family: Insights on whole genome duplication scenarios.

Vertebrate genome evolution remains a hotly debated topic, specifically as regards the number and the timing of putative rounds of whole genome duplication events. In this study, I sought to shed light to this conundrum through assessing the evolutionary history of the oxytocin/vasotocin receptor family. I performed ancestral analyses of the genomic segments containing oxytocin and vasotocin receptors (OTR-VTRs) by mapping them back to the reconstructed ancestral vertebrate/chordate karyotypes reported in five independent studies (Nakatani et al., 2007; Putnam et al., 2008; Smith and Keinath, 2015; Smith et al., 2018; Simakov et al., 2020) and found that two alternative scenarios can account for their evolution: one consistent with one round of whole genome duplication in the common ancestor of lampreys and gnathostomes, followed by segmental duplications in both lineages, and another consistent with two rounds of whole genome duplication, with the first occurring in the gnathostome-lamprey ancestor and the second in the jawed vertebrate ancestor. Combining the data reported here with synteny and phylogeny data reported in our previous study (Theofanopoulou et al., 2021), I put forward that a single round of whole genome duplication scenario is more consistent with the synteny and evolution of chromosomes where OTR-VTRs are encountered, without excluding the possibility of a scenario including two rounds of whole genome duplication. Although the analysis of one gene family is not able to capture the full complexity of vertebrate genome evolution, this study can provide solid insight, since the gene family used here has been meticulously analyzed for its genes' orthologous and paralogous relationships across species using high quality genomes.